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Amendment  No,  1  thereto. 

The  program  was  sponsored  hy  the  Advanced  Research  and  Project  Agency  through  the  Power 
Branch  of  the  Office  of  Naval  Research  under  the  direction  of  ARPA  and  Mr,  J.  A,  Satkowski 
of  ONR. 

The  report  was  written  bv  H.  E.  Wilhelm. 
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I.  INTRODUCTION 


This  report  describes  the  progress  achieved  on  Contract  ONR  No  nr- 4 1 04  (00)  from  1  February 
through  30  April  1065. 

Experimental  effort  during  the  report  period  was  directed  toward  completion  of  the  new  test 
section,  magnet,  and  heater.  In  the  theoretical  portion  of  the  program.  Ohm’s  law  was  derived 
for  nonisothermal  plasma  with  thermal  diffusion.  This  derivation  is  presented  herein. 


1 


II.  EXPERIMENTAL  INVESTIGATIONS 


TTie  magnet  was  operated  with  the  test  section  in  place.  The  highest  magnetic  field  obtainable 
in  the  steady-state  condition  for  extended  time  periods  was  20,000  gauss.  A  calibration  curve 
was  established  showing  magnet  current  versus  magnetic  field.  The  estimated  error  was  less 


than  57- 


Vacuum  checking  of  the  completed  system  was  accomplished.  A  vacuum  of  ion  microns  can  be 
maintained  with  a  mechanical  pump  when  pumping  the  entire  s\stem.  The  system  is  completed 
and  the  first  helium  checkout  runs  are  scheduled  for  the  near  future. 
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III.  THEORETICAL  INVESTIGATIONS 
OHM'S  LAW  FOR  NONJSOTHERMAL  PLASMA  WITH  THERMAL  DIFFUSION 


i 

In  a  plasma,  two  types  of  thermal  diffusion  forms  can  be  distinguished  *  —  those  conditioned 

by  conduction  hoa*  currents  and  those  conditioned  bv  barodiffusion  heat  currents  (radiation  heat 

2 

currents  are  disregarded).  Thermal  diffusion  forces  are  of  importance  for  these  processes 
which  are  based  on  momentum  transfer  between  plasma  components.  In  this  contribution,  the 
transport  of  electrical  current  in  nonisothermai  plasmas  is  treated  for  an  arbitrary  degree  of 
ionization.  It  is  shown,  among  other  things,  that  the  thermal  barodiffusion  effect  (TBDE) 
changes  quantitatively  the  electrical  conductivity  parallel  and  transverse  to  the  magnetic  field 
and  the  Hail  coefficient.  Concerning  the  kinetic  theory  of  the  TBDE,  it  is  demor  utrated  that 
inclusion  of  the  effect  of  intercomponent  momentum  transfer  on  the  relaxation  of  the  individual 
heat  currents  in  the  components  is  indispensable. 

The  following  considerations  apply  strictly  only  to  infrasonic  barodiffusion  velocities: 

7Msr  <vs  -%)2«  kTsr 

This  condition  is  satisfied  for  most  applications — nsr  -  (mg  mr)/(m  +  mr)  -  reduced  mass 
and  Tgr  -  ^sr  F(Ts/ms)  +  (Tr/mr)l  -  reduced  temperature  of  the  s-  and  r-compcnents. 


^Superscripts  refer  to  references  in  Section  IV. 
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THERMAL  DIFFUSION  EFFECT  AND  ELECTRICAL  CONDUCTION 


Consider  a  plasma  composed  rf  electrons  (e),  single  charged  ions  (i),  and  neutrals  (a)  contain- 
ing  an  electromagnetic  field  E,  B.  The  components  s  =  e.  i,  a  are  separately  assumed  to  be 
in  approximate  local  statistical  equilibrium  *  but  relative  to  each  other  in  thermal  nonequilibrium, 
Ts  i  Tr.  The  condition  for  the  latter  is  Eeff  Z  Ep,  where  Eeff  is  the  effective  electrical  field 
acting  on  the  electron  comp  nent  and  Ep  is  the  critical  plasma  field.  Under  these  conditions 
the  transport  of  linear  momentum  of  the  s-component,  interacting  via  friction  forces  ~(vr  -  vg) 
and  thermal  dixfusion  forces  —  [(qr/nr  mr)  -  (q*3 / ng  mg)]  with  the  remaining  components  r  , 
is  described  by:  ’ 


n„  rn 
s  s 


[&+vV] 

E  ' 


vg  +  V  .  Ps  *  ns  es  (E  +  vs  X  B) 


-  ns  ms 


-1 

sr 


r  -  s 


( v  -  v  >  +  a  JE- 

s  V  sr-|^r 


sr 


(JL 

l  ns  ms 


(1) 


The  designations  are  standard.  The  possible  combinations  of  interactions  between  the  charged 
and  neutral  particles  in  plasma  are  treated  as  Coulomb  interactions,  C.I.  ^Dcbye  -  radius  D 


as  maximum  impact  parameter;  transport  cross  section  Qf 


sr 


*/es  erV 

2^kTsr; 


•In  A  \  =  D  ( le  e  J  /3k  T  )”M  and  elastic  sphere  interactions,  S.  I.  interaction  radii 

sr*  “sr  1  s  r1  sr  J  ‘a  9  L 


r^  and  rr;  transport  cross  section  Qsr 

.  7,10 

from  the  13 -moment  approximation, 
diffusion  forces  are 


ir(rg  +  rrT 


]• 


The  coefficients  are  taken  throughout 


According  to  the  latter,  the  coefficients  of  the  thermal 


a  ~  -3/5  . ,  , ,  C.  I. ;  a  =+1/5  - S.  I, 

sr  sr 


The  characteristic  reciprocal  times  determining  the  relaxation  of  the  velocity  fields  are: 
3 


T  ”1  _ 

sr  o 


F  kT  sr  Msr 

if  n  u  m„ 

1  ^sr  s 


n ^  Q g *  •  • ,  C .  I. ,  S.  I. 


(3) 


If  radiation  heat  currents  are  negligible,  the  heat  flux  in  the  s-component  is  determined  by  the 
mechanisms  of  conduction  ~VTg  and  barodiffusion  (Dufour  effect)  -»(vs  -  vr); 
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(A) 


x  Si 

.s  s 


V  T 


,  E 


T  *  .  !LSV  <V.  -  V  >  * 

r  I*  s  sr  sr  °  1  J 

_  5  k 

H  2  ms  Ps  T  s 


The  influence  of  the  magnetic  field  on  the  heat  conduction  process  is  represented  by  the  opera¬ 
tor  (B°  *  B  /  B): 


II  „  X2C  tH*1  B°  (B°.  )  (B°  >-  ), 


n”  =  ( i t2)  ,  X2fi  *<u2r2n  ,  xi^  = «  th 

c~r  g  e  *  o  a  a  a*  ~  . 


(5) 


S  S  S’  s  o  S  S 


Note  that  the  gyration  frequency  is  defined  as  dependent  on  the  sign  of  the  particle,  e " 

*■  23  S 

cg  B/m  ,  The  coefficients  of  the  Dufour  effeet  are: 


Ts  4 

3 _ -9  + - 

.  Tsr  3  ln  A 


-  (>  -  Ml  ■■■  C.I.. 
sr  '  1  s/- 


sr 


(G) 


T  .  T 

16  +  3 _ s  _i8 

T  T 

1  sr  J 


....  S.  I. 


Evaluation  of  the  relaxation  times,  r  g,  for  the  heat  fiuxes  in  the  components  according  to  the 
13-moment  method  yields,  considering  intercomponent  as  well  as  viscous  momentum  transfer;* 


1  4  i  1  ■  ■  i 

*  =  -—  r  1  +  —  r  -1  +  3 


ly 

+ 

aa  ai 


ae+  [Aa] 


+  _  r-  +—  r-.l  +  ] 

e  5  ee  ,0  ea  10  ei  [  ej 


-1  4  -1  13  ,  13 

-  _  T  1  4-  —  f  X  4-  _ *  t  A 

ei 

J  ♦[**]■ 


(7) 


-1  4  -1  10  -l 

T  ,  1  - -  T  1  +  _  T  1  +  3  T 

i  5  u  8  tfl  le 


where  *Sg  and  r  1  arc  given  in  Equation  (3).  The  expressions  Aa,  Ae,  and  A^  stand  for  terms 


*In  Equation  (7),  quantitatively  unimportant  factors,  fgg  =  1  -  (1  -  A  ~g)  (2  lnAgs)  *  1  arising 

from  the  Coulomb  interactions,  are  omitted  for  reasons  of  simplicity. 
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associated  with  higher  approximations  not  considered  here**  Physically,  the  appearance  of  the 
cross  terms  *1  Si  *  (s  i  r)  follows  from  the  interrelation  of  conductive  and  barodiffusion  heat 
fluxes,  Equation  (4). 


In  contrast  to  the  relations  established  in  Equation  (7),  it  was  assumed  in  Reference  7  that  only 

collisions  between  like  particles  are  relevant  for  me  considered  relaxation  (*  "= —  * 

s  5  ss 


l). 


Tliis  approximation,  however,  is  insufficient  as  is  proved  later. 


Assuming  the  Debye-radius,  D,  to  be  small  compared  to  the  dimension  of  the  phisma  system, 
to  which  the  results  shall  apply,  space  charge  effects  can  be  neglected  throughout,  ne  =  n^,  and 
an  ionization  degree  defined  by  «  =  n^K nQ  +  na)  =  ^/(n^  +  na)  .  Equation  (1),  taken  for  s  =  e 
and  s  -  i,  yields  two  independent  equations.  From  these,  one  readily  obtains  alter  elimination 
of  ve  and  by  means  of  the  relations:* 


v*e  =  “v  +  (ve  -  v^)  +  (l  -  "  )  (vj  -  v^) 

(8) 

and 

Vj  =  7  +  (1-x)  (vj  -  va) 

(0) 

which  follow  immediately  from  the  definition  for  the  plasma  velocity,  v*  =  5^ns  mg  v*s  /]T)ns  mg, 
and  of  (ve  -  va)  by  means  of 


0  =  (ve  -  Vj)  -  (ve  -  va)  +  (vj  -  va) 


(10) 


two  equations  for  the  conduction  current  density,  j ,  and  ion  slip  current  density,  "^defined  by: 


j  s  ee  (ve  -  v*)  and  I  =  ^  ci  Ro  -  v^) 

These  equations  are- 

[-  Qe  f  +  Re  (J  •  B6)  B°  +  Se7  X  B°]  +  <rc  j^E+vXB-  E^J 


T  ei 

Tei  +  Tea  - 


Kp  I  +  L_  (I  •  B°)  B°  +  M„  I  X  B° 


(ID 


(12) 


*The  corresponding  equations  used  by  Cowling**  are  applicable  only  for  small  ionization  degrees, 

k«1. 
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-  Qi  j  ' 


t  R  (i  •  B°)  B°  +  S,  i 

1  ■'  r 


X  B«J  +  tr1  [e  +  "v  a  B  -  Ed{  = 


,iL  -  K  1  +  L.  (I  •  3°)  it0  i  M.  I  X  B° 

l  i  i 


In  Equations  (12)  and  (13),  the  diffusion  fields,  caused  by  the  inhomogeneities  in  the  plasma 

1  r  *-♦  ^ 

components,  are  in  nonviscous  approximation  IV  •  Ps  Vp3,  onc^  bemuse  of  me  «  a, 

dv^/dt  "  (-Vp  +1*  X  B)/Xns  gi'.'m  by  (s  =  e,  i): 


r  m 

K-S. 
rn  | 


vp  -  VPs 


**K 


ms  V  T‘l  »  Z*L. 

-  >Lj  sr  sr  i,-t' 

es  ,,  k  sr 

r  X  s 


;«a^rs  bJh-lIz 

niy  nr  m  r 


Further,  <rf,  and  on  designate  the  common  conductivity  coefficients  oi  the  partially  and  f un¬ 
ionized  plasma,  respectively, 

°-o  =  ne  ee2/me  (  Tei  +  7  Jj  <l&> 

*i  =  “i  ei2/mi  Tii  (tf,) 

The  expressions  K  ,  . . .  Sg,  (s  -  e,  i)  arc  given  by,  if  one  makes  the  convention  that  for  s  =  p, 
v~  i  and  for  s  =  i,  e: 


^s  “  ’  ’  +bS^°»-Jl  55  4  es^ 


Ra  =aen"  *b  H 


s  s*  s 


V  y  s 


Ks  2  1  ‘  Ws^s  ’  hsn°^  s  f  V 

Ls  =  dsn^  -hRn^s 


Ms  -  dsn^  -  hgXl  ^  s  t  (!-*>•«  rsn 


S  -  a  XI  +  b  XI"?"  +  <a  (  t  .  ■+  T  ) 

o  e  e  e  i  e  el  ea 


Sj  =  a^f  t  ^ftp  -  *  rie 
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t4llwon 


where  for  an  application,  note  that  some  of  the  coefficients 

*-  I  2I 

if  order  (me/m-  )  and  (nie/m^  &) 


in  Equations  (18) 
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I  *  Elimination  of  the 


Equations  (12)  nnd  (13)  ni*e  readily  solvable  for  the  ion  slip  density, 
latter  and  introduction  of  the  operators  (s  o,  i), 


JP  s  '  Kg  f  us  ^  (B°  •  )  -  Ma  (B°  X  ) 


and  the  abbreviations  (s  :  e,  i) 


us  ■  <KS  Ls  +  mJ)  /(K„  -  1,„) 

y  lia  4  Z 

mi  Tei  +  rea 

lead  to  a  generalized  Ohm's  law  in  the  form: 


Partially  Ionized  Plasma,  0<K  -  1: 


%  {Ke  +  V  +  7  Ql  (Kl  +  LV  j 


-  jRe  (Ke  +  ue)  -  r  r£  (Ki  +  Uj)J  (T-  B°)  n° 

-  J(Ke  sf.  -  Qe  Mc)  +  7  (Kt  st  -  Qt  Mj>]  7*  T?J  (23) 

-  f(Me  Sp  -  Qp  Ue)  +  7  (Mj  Si  -  Qi  U.)  J  (j* X  B°)  X  RO 

Equation  (23)  indicates  that  the  TBDE  modifies,  besides  the  operators,  Y'..,  all  current  terms, 

1 2  r  * 

in  particular  those  assoicated  with  the  Hall  effect  lof  importance  at  free  electron  spiraling— 

i  i-i  «  -l  ^  1  L  u  r 

i.e.  ,  /«e  6>|( Tp.  +  r^)  ~  1  land  the  ion-slip  effect  of  importance  at  free  electron  and  ion 

spiraling — i.e.  ,  l*>  l(*  ■  +  r  nt)  t  •  ^  il  ,  whiie  the  term  ~  (i  .  B°)  B°  is  due  alone  to 

thermal  barodiffusion.  For  explanation,  it  is  referred  to  the  structure  of  the  operator, 
detrrming  the  heat  fluxes  and,  thus,  the  thermal  diffusion  forces  in  magnetic  field.  The  various 
expressions  -- 7 are  not  negligible  for  strong  nonisothermal  plasmas,  Tes  »Tsr,  (s,  r  4  e),  £in 
Equation  (22),  »  ^  =  V<m;  12  mj  (T^/T.J  Qea  /Q  ja  and  ’J 

(n -  Q  .)/(n  Q.)  especially  if  the  ionization  is  appreciable  (ion  slip  effect), 

*  i  3  lo  I  *" 
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In  most  experimental  situations  the  inhomogeneities  in  the  components  are  sufficiently  small 

-*d  — * 

to  allow  regarding  the  diffusion  fields.  Eg,  as  unessential  compared  to  the  applied  fields,  F, 

and  induced  fields  v  X  B.  In  this  quasihomogeneous  approximation,  E“  -  0,  one  can  deduce 

from  Equation  (23)  simple  relations  for  the  current  density  parallel  (II  )  and  transverse  (1)  to 

— *  — »  — *„  — » 

the  magnetic  field  by  splitting  the  field  vectors  in  II  and  A  components,  j  q  =  (j  ■  B  ) 

-  j  =  b°  X  (fx  B°),  etc.  Thus,  one  finds  for  the  current  transport  along  the  B  field  lines 


(24) 


where  the  'parallel  conductivity '*  is  given  by  |^see  also  Equations  (21)  and  (22)J  : 


e(Ke  +  Ue)tT <r.  (K^+  Uj) 


!•  (Qe  -  Re)  (Ke  +  Ue)  +  7  (Qi  -  IT)  (Ki  +  U.) 


(2  5) 


Because  Q  -  R  "  Qe  (B  -  0)  and  K  +  U  «  (K?  +  M?)/K  (B  =  0), 


o-  e  f 

r,{  "  Qc  (B-0) ' 


1  +• 


mp  Tia  K„  (B-0) 


mi  Vs  Ki  (B=0) 


1  (  Tia 

.  mi  v,i*WKi 


R  (B=0)  Q,-  (B-0) 


1-1 


x  (B=0)  Qe  (B  =  0)J 


(2  B) 


Consequently,  o-,f  is  independent  of  the  magnetic  field.  Under  conditions  where  f  =  i,  it  is 
t r  >  as  0  <  Qe  (8=0)  -  1  -  (art  +  bp  +  cp)  <  l.  Note,  however,  that  for  strong  nonisothermy, 
TeS  »Tgr  (s,  r  i  e)  and  appreciable  ionization,  f  §  1.  In  the  same  way,  one  finds  for  the  current 
transport  transverse  to  the  B  field  lines: 

(27) 


A  (E 


v  X  B)  -  0  (E  +  v  X  B)  X  B° 


where  the  transverse  conductivity  is  given  by  see  also  Equations  (21)  and  (22)  : 


<re  Ke  *  ^;K1)X-(,eM(i  +  y«r.Mi)Y 

A - .A - * - 4 - - 


X“  +  Y2 


(28) 


with 


X  5  (Ke  Qe  +  Me  Se)  +  y  (K|  Qi  +  M-  ST 
Y  S  (Ke  Se  -  Me  Qc)  +  7  (K*  Si  -  Qi) 


(29) 
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Thr  quantity  fi  has  the  moaning  of  an  effective  Hall  coefficient  given  by 


Ki  Ke+  y(riKi>Y^ffc  Mc^r«r.  MilX  >(} 
(»0Kei  y  o-j  Kr)  X  -  (>_  Me  17  «T  .  Mo  V 


(30) 


The  complexity  of  the  expressions  for  tr  and  £  is  caused  primarily  by  the  7  terms  discussed 
previously.  Suppressing  the  latter  results  in: 

v  -  (%/qp)  u  +rV\  0*  -  -se/Qc. 

-  I  _  ]  ~  1 

Suppressing  in  these  relations  the  TBDE — Q0— ♦!»  K  — *1,  S  — *  <*»  (r  i  r  )  — lends  to  he 
* r  ■  ”  o  c  e  ei  en 

expressions  of  the  elementary  theory: 


P"*)  ,0  •*-  k 


(  r'1  t  r  ‘V1 
ei  ea 


Comparison  indicates  that<^'  >  o^"  and  p'>P"  for  small  magnetic  fields  with  fle  =”  1.  while 
for  large  magnetic  fields  with £i  <<1,^.1 '  and  fi1  ^  P  '  (reduction  of  the  transverse  TBDE 

by  strong  magnetic  fields). 


From  the  selected  collisional  interactions,  the  elastic  sphere  model  is  much  less  favorable  far 
the  TBDE  (e.g. ,  it  can  be  shown  that  for  predoming  Hug  electron  neutral  interaction,  the  TBDE 
affects  the  conductivity  by  less  than  7.  7%)  than  the  Coulomb  interactions.  For  a  more  quantita¬ 
tive  illustration,  the  electrical  conduction  law.  Equation  (23),  is  now  discussed  for  the  im¬ 
portant  case  of  a  plasma  with  Coulomb  interactions  Tk^.  =  1 ,  Lg  =  =  Us  -  0,  y  -  o"l : 


The  coefficients  <r_,,  Qp,  R^,  and  the  diffusion  field,  Ec  ,  are  obtained  from  the  correspond¬ 
ing  e*  -  cssions  without  asterisk  bv  setting  in  them  terms  and  r_l — interactions  with 

sa  as 

neutrals  a  —  identical  with  zero.  Equation  (31)  exhibits,  besides  the  conductive  thermal  diffu¬ 
sion  effort,  the  TBDE  mid  the  Hall  effect,  where  because  of  =  1,  a  scalar,  the  effective 
electrical  field  acts  directly  on  the  plasma.  It  follows  for  the  current  transport  along  the  B 
field  lines: 
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Allison 


3 


E  -  E 


?*] 


(32) 


where  the  parallel  conductivity  is  independent  of  the  magnetic  field  and  given  by: 
* 


cr  * 


Q*  -  Rf  Qf  (B=o) 


(33) 


Under  consideration  of  Equations  (7)  and  (17)  and  (18),  one  has  explicitly 
2 


-fc) 


« l  T  .  *  T  =  T 

iel  tee  ie 


:  cr*  ~tr*  (13  +  4V2)/4(1  +V2). 


Accordingly,  <r  *  =  1.92  this,  about  twice  as  high  a  conductivity  value  in  comparison  to 
tr*,  agrees  well  with  that  derived  by  Spitzer  by  taking  into  account  the  "electron-electron-in¬ 
teraction.  ”*  One  finds  for  tlie  current  transport  transverse  to  the  B  field  lines: 


cr* 

J  1  “  1 


(E  +  vXB  -  e£)  -  fi  *  (E  +  ~v  X  B  -  E^)x  X  B° 


Tlie  transverse  conductivity  and  the  effective  Hall  coefficient  are  given  by: 


*  /n* 

(r.  t  <re'0 


e 1  <c  R  *  - 

TTpZ’  p 


5L 


>o 


(34) 


(3  3) 


Under  consideration  of  Equations  (7)  and  (17)  and  (18),  one  has  explicitly 


b*  ~  «  1,  T  .  =  T  =  T 

6  \mi  I  1  Cr 


#  •  * 

:  v,  =  tr 
-1  e 


13  t  4>/2 


(1 


4  r*V  1 


e  e 


.  (1  +2*2)f  ,  and  =  +  j  <->L,|  rGi 


1  - 


2  ,  *2v“l 


13  +  4V2 


(1  +  6)  "  T  *«-) 

'  e  c  ■ 


r 


1  + 


90 


- o  (1  t-0)2  r  *2)-  1 

(13  +  4Vr2)2  c  e 


It  is  seen  that  the  TBDE  produces  an  increase  of  the  Hall  coefficient.  This  effect  is,  however, 
the  smaller,  the  larger  the  magnetic  field  intensity  —  e.g.,0  *  =  2.32  j  we  j  Tpj  for  <*»e2  t  ^2«1, 
while  fi*  ~~  |  <up  |  fpj  for  a»G  r^2  »l.  The  transverse  conductivity  is  increased  by  the  TBDE 

for  small  magnetic  fields,  «j2  t  *2  «1,  while  for  large  magnetic  fields,  r  *2  »1,  it  has 

influence — e.g. =,1.92o‘cv  (1  +  «G  re“|)  *  for  a  2  re*2«l,  while  crj*  =  (1  t  «c2  rfc,i‘')  * 


no 


for  w2  r  G2  »1. 
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FURTHER  REMARKS 


It  is  noticeable  that  by  approximating  the  relaxation  times  of  the  heat  fluxes.  Equation  (?),  by 

5 

the  self-collision  terms,  *s  =  —  *_»  *  in  tiie  absence  of  strong  magnetic  fields,  Qp  becomes 
smaller  and  smaller  when  approaching  the  region  of  predominating  Coulomb  intcrations  — 
finally  zero — and  is  negative  throughout  the  region  of  full  ionization.  Completely  analogous 
aspects  are  exhibited  by  Qr  when  the  electron  density  becomes  smaller  and  smaller.  As  Qc 
compares  the  opposing  effects  of  intercomponent  friction  and  thermal  barodiffusion,  this 

approximation  would  lead  to  physically  unacceptable  conductivity  values  (<r— _ L  .' )  and  is, 

consequently,  quantitatively  and  qualitatively  unacceptable. 


As  the  TBDE  and  the  intercomponent  friction  effect  are  opposing  each  other  and  are  for  certain 
collision  models  of  the  same  order  of  magnitude,  their  comparative  consideration  allows  in  a 
simple  way  checking  the  sufficiency  or  insufficiency  of  any  kinetic  approximation  made.  For 
example,  such  comparisons  might  be  useful  in  validating  the  higher  approximations  of  the 
Chapman- Enshog  and  13- moment  methods. 


^Corresponding  to  the  first  approximation  of  Chapman-Fnshog. 
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